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SECTION  I 
INTRODUCTION 


The  rate  of  convergence  of  an  adaptive  filter  algorithm  to  a 
neighborhood  of  the  optimum  filter  is  approximated  in  terms  of  the  eigenvalues 
corresponding  to  the  principal  components  of  some  matrix  closely  associated 
with  the  noise  matrix  (Brennan,  1971).  For  the  frequency-domain  maximum- 
likelihood  adaptive  algorithm,  the  adaptive  update  equation  is: 


.t+1 


I  -  2p 


VV 

V 


Xt(XH)t 


t  th  H 

where  X  is  the  input  data  vector  at  the  t  iteration,  V  is  the  beamsteer 

filter,  A  is  the  conjugate  transpose  of  the  adaptive  filter  vector  at  the  (t+l)1^ 
th 

and  t  iterations,  respectively,  and  fi  is  a  real  scalar  quantity  controlling 

the  adaptation  rate.  The  superscript  H  denotes  conjugate  transpose.  The 

output  of  the  adaptive  filter  at  the  t^  iteration  is  (A^St*.  If  we  substitute 
H  t  H  t 

4>  =  E  (XX  )  for  X  (X  )  in  the  update  equation  for  this  particular  algorithm, 


the  time  constant 


7  is  approximately 
P 


for  energy  lying  on  the 


orthonormalized  eigenvector  of  the  matrix 


vvH\ 

vHv/ 


where  p  <<  1/A  .  Here  A 
P  P 


,  .  th 

is  the  p 


eigenvalue 


of  that  matrix. 
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7  is  the  time  constant  for  the  p  principal  component  in  the  sense  that  the 
P 

portion  of  the  excess  filter  output  RMS 
”\4H4>A  -  (AH«1>A) 

minimum 

th 

associated  with  the  p  principal  component  of  the  matrix 


("*)•(-*) 


is  reduced  by  a  factor  e  =  2.71828  in  approximately  1  / (2^tA^)  iterations. 
Similarly,  the  amplitude 


a  -  (a  ) 

Ip  p  optimum 

of  the  difference  between  the  adaptive  filter  and  the  optimum  maximum- 

th 

likelihood  filter  along  the  p  principal  component  is  likewise  reduced  by  a 

factor  e  in  the  same  period  of  time.  It  is  assumed  that  substitution  of  <i>  for 

t  t  H  j}e 

X  (X  )  does  not  change  the  answer  much. 

Stability  of  the  adaptive  algorithm  requires  that 


P  <  1/A 

max 

so  that  the  time  constant  for  weak  components  may  be  very  great  and  n  can¬ 
not  be  chosen  sufficiently  large  for  practical  effectiveness'  against  very  weak 
components.  The  purpose  of  the  present  study  is  to  investigate  the  ability  of 
modified  adaptive  algorithms  to  increase  the  convergence  rate  on  weak  noise 
components. 

For  the  purpose  of  simplicity,  the  modified  algorithms  in  the 
frequency  domain  will  be  discussed,  although  in  the  following  section  the  mod¬ 
ified  algorithms  are  developed  in  both  the  time  and  frequency  domains.  Let 
there  be  a  sequence  of  data  vectors  X*  at  a  specific  frequency.  These  vectors 

*  This  would  be  the  case  if  successive  data  vectors  were  independent.  In  the 
dependent  data  vector  case  at  hand,  one  might  expect  this  to  be  true  for  /i 
very  small. 


are 


to  be  processed  adaptively  so  as  to  preserve  a  signal  S  =  cV/(V  V), 
where  c  is  a  complex  scalar  quantity  and  V  is  the  conjugate  transpose  of 
the  beamsteer  filter  V  .  One  way  to  formulate  this  algorithm  is 


°1  =  V  X  -  (Fi)X  • 

t  H  t 

where  o  is  the  output  of  the  first  loop  at  time  t  ,  and  (F  )  is  an  adaptive 

*  t  ^  H  t 

filter  applied  to  X  to  predict  the  complex-valued  scalar  quantity  V  X  .  If 

H  t 

the  constraint  (Fj  )  V  =  0  is  imposed  for  all  t  ,  then  it  follows  for  the  input 
Xt  =  S  =  cV/(V^V)  that  the  output  is 

t  cVHV  -  (FjVvj 


vHv 


cVHV 

vHv 


=  c 


H  H 

so  that  the  combined  filter  V  -  Fj  has  a  unit  response  to  the  signal.  It  is 
possible  to  consider  a  second  loop  or  filter 


t  t 

°2  =  °1 


(f^v 


or  in  general 


t  t  ,vFtt 

°k  =  °k.i  -  (Fk'x 


where  the  signal  is  still  preserved  if  the  constraints 


(fJVv  =  0 

K 
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■ 


are  required.  If  the  filters  are  to  operate  independently,  the  optional  addi¬ 
tional  constraints 

(FH)1  F*  =  0  (  j  *  k) 

J  K 

need  to  be  specified.  The  algorithm  with  only  the  constraint  to  preserve  the 

signal  will  be  termed  the  single -constraint  mode  algorithm,  whereas  for  the 

algorithm  with  the  additional  constraints,  the  terminology  multiple -constraint 

H 

mode  will  be  used.  Heuristically,  the  second  loop  allows  the  filters  F  to 

H  ^ 

concentrate  on  the  smaller  components  of  the  noise  because  F^  has  reduced 
the  largest  noise  component.  The  following  section  will  develop  the  algorithm 
mathematically,  and  then  in  the  final  section  some  results  will  be  given  to  in¬ 
dicate  the  behavior  of  the  algorithms  on  synthetic  data. 


1-4 
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SECTION  II 

MATHEMATICAL  DEVELOPMENT  OF  THE  TECHNIQUE 


A.  MULTILOOP  ADAPTIVE  FILTERING  -  FREQUENCY  DOMAIN 

1,  Steepest  Descent  Algorithm  with  Complex  Linear  Constraints 


Let  f  be  a  real-valued  function  of  the  real  variables  Xj,  y^ 

x  ,  . . x  ,  y  .  Suppose  that  these  variables  are  subject  to  the  m  com- 

2  *  2  n  n 

plex  constraint  equations  g^(xj»  Yj»  X2*  Yg*  ••••  XM"  Yr )  =  0  (k  =  1»2 . m). 

To  reduce  the  quantity  f  using  the  steepest -descent  algorithm  while  satisfy¬ 
ing  the  constraint  equations,  the  variables  Xj,  y^  x^,  y^ . y^  can 

be  altered  according  to  the  equation 


where  the  superscripts  t  and  t+1  refer  to  the  t  and  (t  +  1)  iterations, 

respectively,  where  p  is  a  scalar  controlling  the  magnitude  of  the  change  in 

the  vector  P,  where  the  variables  (A  ),  and  (A  ).  are  real-valued  Lagran- 

X  k  y  K 

gian  multipliers  corresponding  to  the  real  and  imaginary  parts  of  the  m  com¬ 
plex  constraint  equations,  and  where  the  vector  P  and  the  operator  V  are 
column  vectors  given  by  the  equations: 


and  V  =  I  d/d  y. 


d/dx 

r 

d/d  y 


The  update  equation  implies  that 


r--;  -g  [<* 


<MRe  8k)  d(Im 

x^k  5xj  +  ^y^k  dx 


yt+1  =  *  .  M  )  *L  +  2 

j  j 


g[“ 


d^e  8i)  dO"i 

A  dVj  +{Ay,k_ d7 


Alternatively, 


#  i  •  ,t+l 

(xj  +  »y.) 

«/  J 


Viy/  -*((%  +  'w) 


I"  d(Re 

S  (Ax)k  ~^7 

k  1  L  i 


.  3(Rc  gk) 


f3«lm  .J  d(Im8,.)l)1 


For  linear  constraint  conditions  of  the  type 


Si 


(Rc  Vk)  -  l  (hr. 


Vj]  '‘j  +  ‘yj> 


■S-. 


(R«  Vj  +  y^Im  VJ, 


k’jl  +  ![yj(R'  Vj  *xJ0m  Vj] 


=  0, 


d(Re  g  )  d(Re  g  ) 

+  1  S7T  =  (Re  Vj  +  1  (Im  Vj 


j 


'j 


and 


j 


'j 


d(Im  g  )  d(Im  g  )  r 

■5^—  *  1  IT  =  +i  L(Re  Vj  +,<Im 


Vj]  ' 


so  that 


,  ,  ,t+l  .  ,  tt 

(*j+*Yj)  *  (*j  +  iYj)  -  P 


i(%  +i^) 


m  r 

E  k 

k=i  L 


tz  >:  "Vk  +  *<Vk]  [(R'vk’i  +  “,mVj 


ii 


where  is  an  arbitrary  complex-valued  column  vector  with  n  complex¬ 
valued  components  (Re  V  )  +  i  (Im  V  )  ( j  =  1 ,  2 . n).  Setting  A  =  (A  ).+ 

k  J  K  j  k  x  k 

1  (A  ).  an<*  (V  )  =  (ReV  ),  +  i  (Im  V  )  ,  the  update  equation  becomes 
yk  kj  kj  kj 


.  ■  ,t+l 

(*j  +  »Yj)  =  (x 


i*V  ’{(£; “■'"■’.I! 
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The  complex -valued  Lagrangian  multipliers  A^  are  determined  from  the 
m  constraint  equations. 

2.  Single  Constraint  Mode 

Let  the  output  o *  at  the  nth  loop  of  a  multiloop  frequency 
domain  adaptive  algorithm  be  given  by  the  formula 


t 

o 

n 


-  (F^Vx1 


where  o*  is  the  output  of  the  (n-l)*^  loop,  F  is  the  n  -loop  adaptive 
filter,  and  X  is  the  input  data  vector.  The  superscript  H  denotes  conjugate 
transpose.  The  superscripts  t  and  t+1  denote  the  vectors  at  iterations  t 
and  t+1,  respectively.  From  this  point  on,  the  superscript  t  will  be  under¬ 
stood  if  the  superscript  t+1  is  not  used.  In  the  first  loop  of  the  adaptive 
filter,  the  zero-th  loop  output  is  defined  to  be  the  beamsteer  output  V  X, 
where  V  is  the  beamsteer  filter. 

The  squared  amplitude  °n°n  the  n*  loop  output  is 


* 

o  o 
n  n 


o 


n-1 


YH 

o  ,  X  F 
n-1  n 


♦  fh  xxhf  . 

n  n 


where  the  superscript  *  denotes  complex  conjugate.  In  the  single-constraint 
mode  of  the  adaptive  algorithm,  this  quantity  is  reduced  subject  to  the  constraint 

vhf  =  0 

n 

using  the  steepest-descent  algorithm. 

-2  Re  (Xo*  .).  +  2Re(XXHF). 

n-1  j  n j 


Since 


J(0n0n> 
»<**  Fn>j 
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Il~  !■«  l_  1  ma 


i 


99MMM W&WiW 


d<°  °J 


n  n 


d(IrnF  ). 
n  J 


-2Im(Xo*  .).  +  2Im(XXHF  ). 

n-1  j  n j 


and 


d(o  o  ) 


n  n 


d(Re  F  ). 
n  J 


d(° 


+  i 


n  n 


d(ImF  ). 
n  j 


=  2  {-X  o*  +  XXHF  ) 

n-1  n  j 


where  j  denotes  the  j1*1  component  of  the  vectors,  the  update  equation  is 


,t+l 

n 


=  F  +2|i(Xo*  .  -  XXHF  -  AV) 
n  n-1  n 


=  F  +  2|i  (X  o  -  AV). 
n  n 


H  t+1 

To  determine  the  Lagrangian  multiplier,  the  relations  V  F  =0  and 

H  t  n  H 

V  F^  -  0  are  used,  and  the  filter  update  equation  is  premultiplied  by  V  : 


0  =  VHFt+1 
n 


VHFl  +  2fi(VHXo*  -  VHV) 

n  n  ' 


=  2|*  (VHX  o*  -  AVHV); 

n 


AVHV  =  VHXo*  , 
n 


A  = 


H  * 
V  Xo 

n 


vHv 


Substitution  of  A  into  the  update  equation  yields 


,t+l 


n 


F  +  2ft 
n 


Xo 


n 


VVHX 

VHV 


<] 
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F  +  2|io 
n  n 


I  - 


vvH 

vHv 


X 


vv 


H 


vHv 


F  +  2|io  X 
n  n 


In  the  last  step,  the  fact  that  V  F  =0  was  used. 

n 

3.  Multiple  Constraint  Mode 

In  the  multiple -constraint  mode,  the  output  o*  at  the  n^ 

n 

loop  of  the  algorithm  is  again 
t  t  H  t  t 

°,  =  °n  1  -  (FJ  X  . 

ii  n-i  n 

but  the  squared  output  °n°n  reduced  subject  to  the  multiple  constraints 


Fk  1  F  =  ° 
k-1  n 


(k  =  1 ,  2,  ,  .  .  ,  n). 


H 

where  F  is  the  beamsteer  filter.  Therefore  the  update  equation  is 
o 


Ft+1  =  F 
n 


+  2p  (Xo  -  F,  A  ,)  .. 

n  n  k-1  k-1 

k=l 


If  the  matrix  M  and  the  column  vector  A  are  defined  by  the  equations 


M  = 


vH 

A 

F» 

o 

A 

1 

1 

• 

and  A  = 

• 

• 

• 

• 

• 

fh, 

A  , 

n- 1 

n~  1 
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the  update  equation  may  be  rewritten 


FtH  =  F  +  2m  (Xo*  -  MHA). 
n  n  n 

To  determine  the  column  vector  A  of  Lagrangian  multipliers, 
the  update  equation  is  premultiplied  by  M*"*"  : 


^  wt+l  ^t+1 
O  =  M  F 

n 

=  Mt+1  F  +  2p(o*Mt+1X  -  Mt+1MHA)  , 
n  n 


where  O  is  a  column  vector  with  n  zeroes; 


A+l  „H 


*  ,  ,t+l 


2  a  M  M  A  =  M  F  +  2p  o  M  X  ; 


a  1  x  *  /AAt+1  aaH\-1  ™t+1  Y 

A  =  — —  (M  M  )  M  F  +  o  (M  M  )  M  X  . 

2  fx  n  n 


Substituting  A  into  the  update  equation, 

t+1  * 

F  =  F  +  2(i  o  X 

n  n  n 


-  MH(Mt+1MH)_1  Mt+1Fn  -  2p  o*  MH(Mt+1MH)_1Mt+1X 

=  |  I  -  (MH)1  ^Mt+1(MH)tJ"1Mt+1J  [V  +  2p(o*)VJ. 


In  the  second  loop, 


_t+l 


I  - 


vv 


H 


vHv 


F*  (f”)'+1 
<f»,'+1f; 


t  *  t  t 

F*  + 


B. 


MULTILOOP  ADAPTIVE  FILTERING  -  TIME  DOMAIN 


1. 


Single  Constraint  Mode 


t  til 

Let  the  output  o  at  the  n  loop  of  a  multiloop  time -domain 
n 


adaptive  algorithm  be 


t  _  t  T  t  t 

n  n- 1  n 

t  til  T  tli 

where  o^  is  the  output  of  the  (n-1)  loop,  Fn  is  the  n  -loop  adaptive 

filter,  and  X  is  the  input  data  vector.  The  superscript  T  denotes  trans¬ 
pose.  In  the  first  loop  of  the  adaptive  filter,  the  zero-th  loop  output  is  defined 

T  T 

to  be  the  beamsteer  output  V  X,  where  V  is  the  beamsteer  filter. 

2  th 

The  squared  output  o^  of  the  n  loop  output  is 


=  o 


n-1 


T  T  T  T 

FXo  ,-o  ,  X  F  +  FXX  F. 

n  n-1  n-1  n  n  n 


In  the  "single -constraint"  mode  of  the  adaptive  algorithm,  the  squared  output 
°n  *s  reduced  subject  to  the  non-r edundant  constraints  expressed  in  the  matrix 
equation 

MF  =  O 
n 

using  the  steepest-descent  algorithm.  Here  M  is  a  matrix  with  as  many  rows 

as  constraint  conditions  and  as  many  columns  as  filter  weights,  F  is  a  column 

n 

vector  containing  the  filter  weights,  and  O  ia  an  all-zero  column  vector  with 
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as  many  elements  as  constraint  conditions.  Several  constraints  were  permitted 
so  that  the  frequency  response  to  a  signal  from  the  desired  look  direction  for  the 
combined  filter 

rri  _£L  rp 

v  -  £ 
could  be  made  white. 

The  update  equation  for  the  adaptive  algorithm  is 


Fn  -  MV 


+  2  A  MF 


=  F  +2// 
n 


T  T 

X  F)  X  -  M  A 


=  F  +  2u  o  X  -  M  A 

n  n 


where  V  denotes  a  column  vector  operator  containing  partial  derivatives  with 
respect  to  the  corresponding  filter  weights  in  F^,  and  where  A  is  a  column 
vector  having  as  elements  Lagrangian  multipliers  corresponding  to  each  of  the 
constraint  conditions.  The  constraint  conditions  MFfl  =  O  and  MFfl  =  O 


are  us 


ed  to  determine  the  column  vector  A  by  premultiplying  the  update 


equation  by  M  : 


O  =  MFt+1  =  MF  +  2u  o  MX  -  MMTA 

n  n  r  |_  n  J 


2  u  o  MX  -  MM  A 

’  rt 


MM  A  =  o  MX  ; 

n 
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A  =  o  (MM1)'1  MX. 
n 


Ft+1  =  F_  +  2m„  [  X  -  MT  (MM1)'1  Mx] 


Substitution  into  the  update  equation  gives 

a 

=  Fn  +  2*.  On  [  I  -  MT  (MM1)'1  m]  X 

=  -  MT(MMT)-1  M  ]  [  Fn  +  2m  on  x]  . 

2.  Multiple  Constraint  Mode 

The  output  o*  at  the  n^  loop  of  the  algorithm  is  the  same 
as  for  the  "single-constraint"  mode,  but  the  constraints 


F,  F  =  0 
k  n 


(k  =  1 ,  2,  .  .  ,  ,  n-1) 


supplement  those  expressed  in  the  matrix  equation 


M  F  =  O  , 
on  c 

where  M^  is  a  time -invariant  matrix  used  to  express  c  linear  constraint 

conditions  and  O  is  a  column  vector  with  c  zeroes, 
c 


The  update  equation  is 


^t+1 
F  =  F 
n 


-  mV  To2  +  2  ATM  F  +  y.  2  A,  F,TF  1 

n  L  n  o  o  n  k  k  nj 


n-1 


=  F  +  2m  (o  X  -  M  A 

n  n  o  o  k  k 

k=l 


23  Fk^k^  * 


where  A^  is  a  column  vector  containing  Lagrangian  multipliers  correspond¬ 
ing  to  the  c  constraint  conditions  expressed  in  the  matrix  equation  M  F^  =  O  , 
and  where  the  values  A  are  scalar  Lagrangian  multipliers  corresponding 

1C 

to  the  new  constraint  conditions. 
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-  •IflfegfeAo. . 


H&s** 


:i 


If  a  new  matrix  M  and  column  vector  A  are  defined  by  the  equation 


the  update  equation  may  be  rewritten 

t+1  T 

F  =  F  +  2  (o  X  -  M  A  )  . 

n  n  n 

M  now  has  c  +  n-1  rows,  and  A  now  has  c  +  n-1  elements. 

To  find  the  column  vector  A  ,  the  update  equation  is  pre- 
multiplied  by  Mt+S 


O 


c+n-1 


..t+l^t+i 

M  F 

n 

Mt+1F  +  2/i  (o  Mt+1X  -  Mt+1MTA)  . 
n  n 


O^n  j  is  a  column  vector  with  c+n-1  zeroes.  The  superscript  t+1  is 

attached  to  the  matrix  M  because  it  contains  the  time  varying  filters  F' 
T 

F^  j.  Rearrangement  of  the  preceding  equation  gives  the  formula 


■y  xit+lwT.  t+l_  ,  x.t+l 

2 ii  M  M  A  =  M  F  +  2/io  M  X 

n  n 


and  hence  the  result 


T  T  -1 
M  (M  M  ) 
o  o  o 


1 


1  ,_T  t+1  t 
|  (F,  )  F 


1 


M 


(F^)t+1 


Fl  (FT)t+1 

T  T-l  l  '  1  # 

=  M  (M  M  )  M  +  - - 

o  o  o  o  ,„T.t+l  _t 


(F j  ) 


1 


Hence  the  second-loop  update  equation  is 


,t+l 


T  T-l 
I  -  M  (M  M  )  M 
o  o  o  o 


f\  (F^)t+1 
(F^)t+1  f\ 


[f2  + 


c. 


DATA  GENERATION 


An  array  of  elements  located  in  a  three  dimensional  space  is 
defined  by  a  set  of  coordinates,  (X^,  Y^,  Z^)  where  k  =  1,  2,  ,,,,  C,  and 
C  is  the  number  of  elements  in  the  array.  Signal  vectors  in  the  frequency 
domain  as.d  at  frequency  f,  S^,  for  the  array  are  defined  for  each  of  the  spec* 
ified  "lock  directions"  by  the  equation: 


S.  = 
J 


-i27T  f  A  tj  (1) 


~i2.1T  f  A  t.  (2) 
I  J 

l 

i 

-i2  77  f  A  t.(C) 


where  j  references  a  specific  "look  direction".  The  j**1  "look  direction"  is 

defined  by  a  set  of  angles  AZ.  (Azimuth)  and  EL.  (Elevation)  as  illustrated 

J  1 

below. 
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The  A  t.(k)  is  the  relative  time  delay  of  the  signal  at  array 
element  k  and  is  dependent  upon  the  "look  direction",  j,  and  the  propagation 
velocity  of  the  signal,  VEL. 


Consequently 


S. 

J 


S, 

J 


can  also  be  defined  as 


U-  R(l) 
VEL, 

U-R(2) 

VEL 


U « R(C) 
"  VEL 


where  R  (k)  is  a  vector  in  the  array  coordinate  space  defining 
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the  location  of  element  k  relative  to  some  reference  element,  and  U  is  s  unit 
vector  from  the  reference  element  in  the  direction  of  the  "look  direction".  There¬ 
fore,  the  components  of  U  are: 


Ux  =  Cos  (EL)  Sin  (AZ) 
Uy  =  Cos  (EL)  Cos  (AZ) 
=  Sin  (EL) 


A  composite  cross-power  matrix,  O  ,  is  generated  from  p  signal 
vectors  plus  random  noise  according  to  the  following  formula: 


r 

n*  X  s(sf  ♦  ei 

j=l  J  '  J 


where  a.  and  $  are  the  relative  weights  to  be  associated  with  the  U>lgnal  vectors 
J 

and  random  noise  component,  respectively. 

If  we  let  x  be  v  complex  random  vector  such  that 


P 

_  V 


4—  (|  V  4 W 

j-1  •  J 


where  (.  is  a  random  number  with  normal  distribution,  zero  mean  and  variance 
a  ,  ,  and  where 


w  is  a  complex  random  vector  such  that 


=  1  ,  then  by  choosing 


11-16 


it  follows  that 


\H 

c.S.  +  6wJ 

} 3  I 


p  — 


=  £ 


J=1 


Zc  r-H 

t  s  r 
J  i  J 


+  b '■ 


ww 


H 


■  §  aJsJ5"  ♦  61 


Thus,  we  can  use  equation  (2)  to  generate  a  set  of  data  vectors  which  has  the 
crosspower  matrix  f)  in  equation  (1 ). 
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D. 


EVALUATION  AND  DISPLAY 


We  try  to  evaluate  the  performance  of  the  multiloop  adaptive 
filtering  technique  by  comparing  its  mean  square  output  (MSO)  at  time  t  for 
each  of  the  signal  vectors,  white  noise  and  composite  energy  with  those  of  the 
optimum  maximum-likelihood  filter. 

The  formula  we  used  to  calculate  the  MSO  is  as  follows: 


Optimum  Filte  r 


Adaptive  Filter  Loop  K 


White  Noise 


Signal  Vector  j 


J 


H  F 

0(V 

o  o 

2 

F  S. 

a. 

°  J 

J 

K  K 

V  „t  ...  V 


i=l 


i=l 


(V 


-  t  1'')HS. 

1  1  J 


Composite  Energy 


fh  n  F 
o  o 


(V  -  Z  f‘  )H  n  (V.  f  F* 
i=l  1  i=l  1 


where  Fq  is  the  maximum-likelihood  filter  and 


F  = 
o 

where  Z  =  fl"*  Sj 


and  Sj  is  the  signal  vector  to  be  preserved. 

We  take  the  ratio  of  optimum  filter  MSO  to  adaptive  filter  loop  k 
MSO  and  plot  it  as  a  function  of  t  (number  of  adaptive  iterations)  for  each  of  the 
signal  vectors,  white  noise  and  composite  energy,  respectively. 
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SECTION  III 

PROCESSING  RESULTS  AND  CONCLUSIONS 


A  set  of  synthetic  data  was  generated  for  an  equally-spaced  line 
array  of  five  elements.  Three  coherent  components  were  added  to  synthetic 
random  noise.  The  first  of  the  three  coherent  components,  S^,  was  a  signal 
from  the  look  direction.  The  second  and  third  coherent  components,  S  and 

M 

S^»  were  interfering  noise  components.  Different  weightings  were  assigned 
to  these  vectors  (as  shown  in  Figure  III- 1 ).  The  weighting  =  0  for  the 
signal  component  S^  indicates  that  the  signal  component  was  not  considered 
to  be  part  of  the  noise  output.  Component  3  represented  a  weak  noise  com¬ 
ponent  in  the  presence  of  a  strong  noise  component  (component  2).  It  was  ex¬ 
pected  that  multiloop  adaptive  filtering  would  increase  the  sensitivity  of  adap¬ 
tive  processing  to  the  weak  noise  component. 

Several  runs  with  two -loop  adaptive  filters  were  made  using 
500  synthetically-generated  data  vectors  as  input  to  a  like  number  of  iterations 
of  the  adaptive -filter  algorithm.  Using  the  frequency-domain  algorithms  de¬ 
rived  in  the  previous  section,  the  adaptive -filter  algorithm  was  run  both  in  the 
single -constraint  mode  and  the  multiple -constraint  mode.  Two  different  sets 
of  initial  filters  were  tried.  The  first  was  an  optimum  filter  designed  for  the 
data  as  specified,  except  that  component  3  was  omitted.  The  second  was  a 
simple  beamsteer  filter. 


The  curves  plotted  are  the  ratios 


Fn  n  F 

O  CO 


A  "  ‘  A  ’ 


where  is  the  crosspower  spectrum  matrix  for  each  of  the  noise  compon- 

H  H 

ents,  Fq  is  the  optimum  filter  and  F^  is  the  composite  adaptive  filter 
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VELOCITY  =  4 
FREQUENCY  = 
WHITE  NOISE  3 


Figure  III- 1 .  Data  Model 


t  mm 


vector  corresponding  to  the  iteration  for  which  the  power  ratios  are  plotted. 
Each  curve,  therefore,  measures  the  performance  of  the  adaptive  filter  rela¬ 
tive  to  the  optimum  filter  for  each  noise  component  or  all  noise  components 
combined. 

Figures  111-2  and  III- 3  show  the  single-constraint-mode  results 
for  the  one-loop  and  two-loop  algorithms  through  500  iterations.  The  starting 
filter  was  an  optimum  filter  designed  for  the  specified  noise  field  buth  with 
component  3  omitted.  Convergence  rates  were  =  0.  00015  in  the  first  loop 
and  n  =  0*  0003  in  the  second  loop.  The  two-loop  algorithm  concentrated 
much  more  effort  on  the  weak  component  and  actually  suppressed  this 
energy  more  than  the  optimum  filter.  The  results  of  Figure  III- 3  reflect  the 
single-constraint  mode  of  processing,  where  the  second-loop  filter  is  con¬ 
strained  only  to  pass  the  signal. 

In  the  top  half  of  Figure  III -4  are  presented  the  numerical  values 
of  the  optimum  filter  set.  Differences  between  the  optimum  filter  set  and  the 
one -loop  or  two-loop  filter  sets  after  500  iterations  are  shown  in  the  bottom 
half.  For  purposes  of  comparison,  the  differences  between  the  optimum  filter 
weights  and  the  initial  filter  weights  are  also  shown  in  the  bottom  half  of  the 
figure.  The  circles  representing  the  two-loop  system  lie  much  closer  to  zero 
than  the  triangles  corresponding  to  the  one-loop  system.  Thus  the  two-loop 
system  appears  to  have  converged  significantly  faster  toward  the  optimum 
filter. 

A  similar  set  of  results  is  shown  in  Figures  III-5  and  III-6.  In 
this  case,  the  convergence  rates  were  increased  and  the  beamsteer  filter  set 
was  used  as  the  initial  filter  se',.  Again,  the  two-loop  algorithm  appears  to 
have  converged  more  rapidly  toward  the  optimum  filter  and  to  have  responded 
more  rapidly  to  the  weaker  noise  component. 

Figure  III-7  gives  the  multiple -constraint  results  for  the  two- 
loop  algorithm  after  500  iterations}.  It  can  be  compared  with  Figure  III- 3  to 

III- 3 
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Figure  III- 5.  Single  Constraint  Mode  -  One  Loop  Starting  with 
Beamsteer  Filter  ( /u,  =  .001) 
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Figure  III-6.  Multiple  Constraint  Mode  -  Two  Loops  Starting  with 
Beamsteer  Filter  (  =  •  001,  =  •  0015) 
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determine  the  effect  of  the  constraint  requiring  orthogonality  between  in 
the  first  loop  and  in  the  second  loop.  Both  figures  depict  results  obtain¬ 

ed  from  an  initial  filter  designed  to  minimize  the  specified  noise  field  (with 
the  weaker  coherent  noise  component  removed).  Convergence  rates  in  both 
cases  were  =  0.00015  in  the  first  loop  and  M  =  0.0003  in  the  second 
loop.  The  orthogonality  constraint  appears  to  have  stabilized  the  convergence 
process.  The  noise  output  power  for  the  individual  components  did  not  fluctu¬ 
ate  nearly  as  much  in  the  multiple -constraint  mode.  Furthermore,  the  noise 
output  power  for  the  various  noise  components  after  500  iterations  wss  much 
closer  to  that  which  would  have  been  obtained  from  an  optimum  system. 

It  is  interesting  that  the  various  adaptive  algorithms  perform 
better  than  the  optimum  filter  on  the  discrete  coherent  noise  components  at 
the  expense  of  poorer  performance  on  the  white  noise  and  composite  noise. 

The  chief  objective  of  this  study  was  to  determine  whether 
multiloop  adaptive  filtering  could  reduce  the  differences  in  the  rate  of  conver¬ 
gence  associated  with  unequal  eigenvalues  in  the  noise  matrix.  The  simulation 
results  have  demonstrated  that  multiloop  adaptive  filtering  does  indeed  have 
this  capability.  This  capability  is  of  potential  use  in  speeding  the  rate  of  con¬ 
vergence  when  the  noise  statistics  vary  slowly  with  time.  In  the  event  of  such 
quasi- stationary  noise  statistics,  the  single-loop  adaptive  algorithm  would 
otherwise  tend  to  hinder  adjustment  to  the  changing  noise  statistics  because  of 
unequal  eigenvalues  in  the  noise  matrix. 
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